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Abstract. The unidentified TeV source in Cygnus is now confirmed by follow-up observations from 2002 with the 
HEGRA stereoscopic system of Cherenkov Telescopes. Using all data (1999 to 2002) we confirm this new source 
as steady in fiux over the four years of data taking, extended with radius 6.2' (±1.2gtat ±0.9sya) and exhibiting 
a hard spectrum with photon index —1.9. It is located in the direction of the dense OB stellar association, 
Cygnus 0B2. Its integral flux above energies E > 1 TeV amounts to ~5% of the Crab assuming a Gaussian 
profile for the intrinsic source morphology. There is no obvious counterpart at radio, optical nor X-ray energies, 
leaving TeV J2032-I-4130 presently unidentified. Observational parameters of this source are updated here and 
some astrophysical discussion is provided. Also included are upper limits for a number of other interesting sources 
in the FoV, including the famous microquasar Cygnus X-3. 

Key words. Gamma rays; observations - Stars: early-type - Galaxy: open clusters and associations: individual: 
Cygnus OB2, Cygnus X-3 



1. Introduction 

The reasonably large fields of view (FoV, eg. FWHM 
> 3°) achieved by ground-based 7-ray telescopes per- 
mits survey-type observations using a few or even singly 
pointed observations. Such potential was realised with the 
serendipitous discovery of a TeV source in the Cygnus 
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region. Analysis of archival data (~121 h from 1999 to 
2001) of the HEGRA system of Imaging Atmospheric 
Cherenkov Telescopes (HEGRA lACT-System) revealed 
convincing evidence for an apparently steady, spectrally 
hard (—1.9 differential photon index) and possibly ex- 
tended TeV source (Aharonian et al. I2002|l . Follow-up 
observations using the HEGRA lACT-System were per- 
formed during its final season of operation (2002) for a 
total of ^158 hours. Analysis of these data again re- 
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veal the presence of this source, thus confirming its ex- 
istence with the HEGRA lACT-System. Earher, observa- 
tions (1991) with the HEGRA scintiUator array (Merck 
I1993|l revealed a multi-TeV excess, positionally consistent 
with Cygnus X-3. Analysis with improved direction recon- 
struction (Krawczynski' 19951 revealed this excess (-1-4.3(7 
pre-trial) as centered roughly 0.5° north of Cygnus X- 
3. Interestingly, the Crimean group (using the Cherenkov 
imaging technique) reported a significant excess (^ +6.0a 
pre-trial) ~ 0.7° north of Cygnus X-3 (Neshpor et al. 
I1995|l . and recently, the Whipple collaboration also re- 
ported an excess at the HEGRA position {+3.3a) in their 
archival data of 1989/1990 (Lang et al. lMnji . 

We summarise here in some detail our analysis and 
numerical results for TeV J2032-I-4130, and give a brief 
astrophysical interpretation. Given the large FoV of obser- 
vations, deep exposures were also obtained for a number 
of other interesting sources. Upper limits for these sources 
are given. 

2. Data analysis & results 

The HEGRA lACT-System, de-commissioned in 
September 2002, consisted of five identical Cherenkov 
telescopes (each with 8.5 m^ mirror area) on the Canary 
Island of La Palma (2200 m a.s.l.). Employing the 
stereoscopic technique, this system achieved an angular 
resolution < 0.1° and energy resolution < 15% for 7-rays 
on an event-by-event basis over the 0.5 TeV to >50 TeV 
regime. Detailed technical descriptions of the HEGRA 
lACT-System and performance can be found in Piihlhofer 
et al. ifMTSji . 

Gamma-ray-like events are preferentially selected 
against those of the dominant isotropic background of 
cosmic-rays (CR). This is achieved with cuts on 6, 
the angular distance between the reconstructed and as- 
sumed arrival directions, and also the mean- scaled-width 
tv (Aharonian et al. I2000a|l . which is a measure of an im- 
age's conformity to a 7-ray-like shape. Event arrival direc- 
tions are reconstructed using the algorithms described in 
Hofmann et al. H1999|l . We a priori selected algorithm '3' 
for final analysis, but also employed the other available 
algorithms to check consistency of the signal. So-called 
tight cuts were implemented, namely: 9 < 0.12°, Hj < 1.1, 
which are optimal for point-like sources in a background- 
dominated scenario. We also required a minimum ntci > 3 
telescope images per event for the 9 and w calculation 
(Aharonian et al. 12002(1 . An estimate of the CR back- 
ground surviving cuts is made using both the template 
(see Rowell 12003 1 and displaced background models for 
consistency checks on the source excess. The displacement 
background model employs different regions in the FoV for 
background estimation using events from ui < 1.1. 

In Tables □ and H 

we summarise details of the TeV 
source which includes the excess significance, the source 
extension agrc, both of which are calculated at the excess 
centre of gravity (CoG), and also the energy spectrum and 



flux. Some results are split according to data subsets, 1999 
to 2001 (dataset §1), and 2002 (dataset §2). 

The CoG and source extension were estimated by 
fitting a 2D Gaussian convolved with the system point 
spread function (PSF) to a histogram of 7-ray-like {w < 
1.1) events binned over a 1° x 1° FoV. The PSF is esti- 
mated from Crab data, giving a value which agrees with 
Monto Carlo simulations of a point source (Aharonian et 
al. I2000bl I2()()4|l . A Gaussian profile suitably describes 
the excess morphology (discussed later). Here, we used 
higher quality events selected according to the estimated 
error in reconstructed direction e (see eg. Hofmann et al. 
I2000|l . This cut on e reduced somewhat the systematic 
differences in CoG obtained from all three available algo- 
rithms '1', '2' and '3'. Simulations and empirical results 
from other point sources of the HEGRA lACT-System 
archive have shown that consistent sensitivities from all 
algorithms are expected for a minimum number of images 
per event ntei > 3. In fact the signal excess significance 
from all data using each algorithm are consistent to within 
Icr, and the CoGs and source extensions agree to within a 
2(7 level. Taking the final CoG and source extension from 
algorithm '3', the respective systematic errors were taken 
using results from the other algorithms. The excess signif- 
icance from combined data exceeds 7(7 (-1-7. Itr), and the 
source appears extended at the > 4(7 level when, conserva- 
tively, subtracting the systematic error in extension. The 
source is termed TeV J2032+4130 based on its CoG. 
The post-trial significance is ~ 6.1a if one assumes the 
trial factors (1100) accrued from the discovery dataset §1 
(Aharonian et al. 12002(1 . The radial source extension is 
shown in Fig. ^ and compared against a number of pos- 
sible source morphologies. We tested three different types 
of intrinsic source morphology: 

— Disc: The source emits at a constant brightness out 
to a radius 0.13° (with zero brightness outside). This 
profile resembles that of a disc. 

— Volume: The source emits at all positions inside a 
sphere of radius 0.13°. 

— Surface: The source emits at all positions within a ra- 
dial band between 0.117° and 0.13°. 

In each case the volume-integrated radial brightness is 
calculated after convolution with the instrument PSF, 
and compared with the measured radial profile of 
TeV J2032+4130. In all cases similar reduced-x^ ~ 1-0 
were obtained, thereby not permitting discrimination be- 
tween these source morphologies. 

Fig. 121 depicts the 2D skymap of excess significance over 
the FoV, using the template model as a CR background es- 
timate. TeV J2032-I-4130 is clearly positioned at the edge 
of the error circle of 3EG J2033-I-4118, and within the 
core circle of the extremely dense OB stellar association 
Cygnus 0B2 fKnodlseder lM)n|l . 

The energy spectrum determination followed the 
method described in Aharonian et al. 1(1999(1 using tight 
cuts (on and id) plus an additional cut on the recon- 
structed air-shower core distance of the event core<200m. 
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Fig. 1. (Top) Events (solid points with error bars) vs. distance from the CoG squared 6^ compared with a background 
estimate from the template model (dashed-dotted line). A convolved radial Gaussian fit F = Ped + P2 exp(— 6'^/ {Pi + 

=0.104° ± 0.020° the intrinsic source size. The PSF width 



a 



p^)) is indicated by the solid line with Pi = cTgrc 
is crpt=0.070° (dashed). The pedestal Ped is fitted separately. (Bottom) Comparisons with various intrinsic source 
morphologies convolved with the PSF (see text). 



Reconstruction of event energies employed the more ad- 
vanced method of Hofmann et al. (2000) which makes use 
of the height of shower maximum to improve the core 
distance determination, and hence improve energy reso- 
lution to < 15%. Both datasets §1 and §2 yielded con- 
sistent power-law fits with a hard photon index. For all 
data, a pure power law explains well the energy spectrum, 
showing no indication for a cut-off when fitting also a com- 
bined power law -I- exponential cutoff term exTp{—E/Ec). 
We estimated, nevertheless, lower limits (99% c.l.) to the 
cut-off energy for a range of (fixed) power law indices. 
Cut-off lower limits of Ec ~ 3.6, 4.2 and 4.6 TeV resuff 
when fixing the power index at values 7 =1.7, 1.9, and 
2.2, respectively. 

In estimating all fiux values we have assumed a 
Gaussian source profile according to the estimated source 
size since the non-pointlike nature of the source is con- 
firmed. Note that previously published integral fluxes for 
combined data (Rowell et al. I200HI Horns et al. I2004|l as- 
sumed a pointlike source. It is also apparent that the event 



rate for dataset §2 is ~80% that of dataset §1, and that 
the integral fluxes (for E > 1 TeV) derived differ by about 
40%. However the statistical errors on the integral flux, 
which are dominated by contributions near 1 TeV, sug- 
gests this difference is not significant (r-^ la). 

We conclude therefore that TeV J2032+4130 exhib- 
ited a constant fiux from 1999 to 2002. Note that the in- 
tegral (and differential) fluxes are corrected for changes 
in the lACT-System response over time (see Piihlhofer et 
al. 120031 for details on the system performance over time). 
In this case, corrections up to the individual run level ac- 
cording to the CR background rate (a reliable, relative 
measure of the detector -|-atmospheric transmission) have 
been applied. 

We also include in Table |21 99% upper limits for a 
number of other interesting sources in the FoV. These in- 
clude the two GeV sources, their possible associated X-ray 
counterparts (as indicated by Roberts et al. 1200 l|l . and 
also Cygnus X-3. 
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Fig. 2. Skymap of correlated event excess significance (cr) from all HEGRA lACT-System data (3.0° x 3.0° FoV) 
centred on TeV J2032+4130. Nearby objects are indicated (EGRET sources with 95% contours). The TeV source 
centre of gravity (CoG) with statistical errors, and intrinsic size (std. dev. of a 2D Gaussian, asrc) are indicated by 
the white cross and white circle, respectively. 



3. Discussion and Conclusion 

Possible origins of TeV J2032+4130 have been considered 
in Aharonian et al. ifMT^ . Butt et al. H2003|l . Mukherjee 
et al. p003|) and Bednarek H2003II . One interpretation in- 
volves association with the stellar winds of member stars 
in Cygnus 0B2, individually or collectively, which provide 
conditions conducive to strong and stable shock formation 
for particle acceleration. Another scenario involves parti- 
cle acceleration at a termination shock, which is expected 
at the boundary where a relativistic jet meets the inter- 
stellar medium. TeV J2032-t-4130 in fact aligns well within 
the northern error cone of the bi-lobal jet of Cygnus X-3 
discussed by Marti et al. H2000I nUUT\i . The existence of 
TeV emission clearly suggests that particles are acceler- 
ated to at least multi-TeV energies. Taken at face value 
the different flux levels (a maximum factor ^20 differ- 
ence) claimed by the Crimean (1.7 Crab E > 1 TeV), 
HEGRA (0.05 Crab E > 1 TeV) and Whipple (0.12 Crab 
E > 0.6 TeV) groups over a period of a decade, would sug- 
gest episodic emission from TeV J2032-I-4130. Explaining 



also the extended nature of TeV J2032-I-4130 as seen by 
HEGRA requires consideration of issues such as the par- 
ticle acceleration site and its distance to that of the TeV 
7-ray emission, particle diffusion and source age. Such is- 
sues would also not in general rule out an extended and 
episodic source. Moreover the TeV J2032-I-4130 emission 
at any given time may be a superposition of more than 
one component (eg. variable compact in addition to weak, 
steady, extended emission). 

For illustrative purposes we have matched the spec- 
tral energy distribution of TeV J2032-I-4130 with simplis- 
tic leptonic and hadronic models (Fig.|3Jl. We assume that 
the TeV emission arises from a single pure population of 
either non-thermal hadronic or electronic parent particles. 
We do not consider here the conditions under which par- 
ticles are accelerated or how they lose energy. Under the 
hadronic scenario the 7r°-decay prediction matches well 
the TeV flux using a parent proton power law spectrum of 
index —2.0 with a sharp limit at 100 TeV. The neighbour- 
ing EGRET source 3EG J2033+4118 (possibly not related 
to TeV J2032-I-4130) should be considered here as up- 
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Fig. 3. Spectrum of TeV J2032+4130 (this work - HEGRA) compared with purely hadronic (Protons E< 100 TeV) 
and leptonic (Electrons E< 40 TeV) models. Upper limits, constraining the synchrotron emission (leptonic models), 
are from the VLA and Chandra (Butt et al. 2003) and ASCA (Aharonian et al. I2002|l . In the model a minimum 
electron energy ^jmin ^ 10* is chosen to meet the VLA upper limit. EGRET data points are from the 3rd EGRET 
catalogue (Hartman et al. I1999|l . 



per hmits on the potential GeV flux of TeV J2032+4130. 
Associated synchrotron X-ray emission would also be ex- 
pected from tertiary electrons (tt^ ... fx^ . . . — » e^), 
not modeled at present, which represent an absolute lower 
limit on any synchrotron emission visible assuming a pure 
electronic scenario. TeV data are matched well by an 
inverse- Compton spectrum (from electrons up-scattering 
the cosmic microwave background) arising from an elec- 
tron spectrum with power law index ~ —2.0 and a sharp 
limit at 40 TeV. The predicted synchrotron emission then 
follows as a function of local magnetic field, constrained by 
the available upper limits at radio and X-ray energies. The 
most conservative synchrotron prediction arises for Bq = 
3/xG, the lowest such B-field expected in the Galactic disk. 
In fact, much higher B-fields (i?o > lO/zG) are generally 
expected in such regions containing young/massive stars 
with high mass losses and colliding winds (e.g. Eichler 
et al. 11993(1 . X-ray results from ASCA (Aharonian et al. 
12002(1 provide constraining upper limits, as do results from 
Chandra (Butt et al. p003 ) . Deeper observations by XMM 
and Chandra will no doubt provide further constraints 
on the leptonic component. Future 7-ray observations by 
H.E.S.S., VERITAS and also MAGIC-II will be vital in 



better determining the spectrum over the ^50 GeV to 
>10 TeV regime. Knowledge of the spectral behaviour 
for E >\Q TeV will convey important information on the 
maximum particle energies and their type. Such informa- 
tion can come from low elevation observations by H.E.S.S. 
for which very large collecting areas are achieved at TeV 
energies. Energy-resolved morphology studies can also be 
performed, allowing conclusions on the diffusion proper- 
ties of the accelerated particles. Overall, TeV J2032-I-4130 
is the only TeV source so far without an obvious multi- 
wavelength counterpart, and is likely the first-discovered 
galactic TeV source which is extended in nature. 

The serendipitous detection of such a weak 0.05 
Crab), marginally extended source with a hard spectrum 
over a long observation time illustrates the power of the 
stereoscopic technique as was employed by the HEGRA 
lACT-System. With a sensitivity at least a factor of 5 bet- 
ter, the next generation instruments will find such sources 
detectable in under 10 hours, or even less for sources with 
steeper energy spectra. 

Finally, we also obtained upper limits (for a steady 
TeV flux) from a number of other source positions in 
the FoV (Fig. Table Ej), including the famous micro- 
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Table 1. Numerical summary for TeV J2032+4130. (a) 
Centre of Gravity (CoG) and extension ctsic (std. dev. of 
a 2D Gaussian); (b) Event summary. The values s and b 
are event numbers for the 7-ray-like and background (from 
the Template and Displaced models, see text) respectively, 
and s — a6 is the excess using a normalisation a. S denotes 
the excess significance using Eq. 17 of Li & Ma <|1983|l : (c) 
Integral events after spectral cuts. 



(a) CoG & Extension (e < 0.12°) 

— All Data (278.3 h) — 
RA Q2000: 20'''' 31™ 57.0' ±6Xt^t ±13.7^y 
Dec 52000: 41° 29' 56.8" il.l^^t^t il.O^^y, 
asra 6.2' ±1.2^tat ±0.9' 



(b) Tight cuts: 6 < 0.12°, w < 1.1, ntei > 3 



Background 








Model 


s b as 


— ah 


S 


— 1999 to 2001 Dataset §1 (120.5 h) — 




Template 


529 2432 0.168 


123 


+5.3 


Displaced 


529 6982 0.059 


119 


+5.4 




2002 Dataset §2 (157.8 h) 






Template 


716 3494 0.168 


129 


+4.8 


Displaced 


716 8510 0.070 


125 


+4.8 




— All Data (278.3 h) — 






Template 


1245 5926 0.168 


252 


+7.1 


Displaced 


1245 15492 0.065 


243 


+7.1 



(c) Spectral Cuts: Tight Cuts + core< 200m 

Energy estimation method: See Hofmann et al. pOHHIl 
Background 

Model s b a s — ab S 

— 1999 to 2001 Dataset §1(120.5 h) — 
Displaced 421 2120 0.143 118 +6.0 

— 2002 Dataset §2(157.8 h) — 
Displaced 552 2999 0.143 124 +5.3 

— All Data (278.3 h) — 
Displaced 973 5119 0.143 242 +7.9 



Table 2. Numerical summary for TeV J2032+4130 con- 
tinued., (d) differential fluxes and events after tight spec- 
tral cuts; (e) Fitted power law; (f) Integral flux. 



(d) Spectral Cuts: Differential points 



Energy 


Flux (E)" 


Flux"" 


s 


b 


S" 


£;(TeV) 




Error (£;) 








— 


1999 to 2001 Dataset §1(120.5 h) — 




1.05 


1.70 


1.24 


130 


678 


+3.0 


1.82 


1.94 


0.91 


107 


531 


+3.1 


3.16 


0.37 


0.27 


57 


323 


+ 1.4 


5.50 


0.30 


0.11 


33 


117 


+3.2 


9.55 


0.09 


0.04 


13 


34 


+2.8 




— 2002 Dataset §2(157.8 h) 






1.05 


19.90 


16.57 


153 


885 


+2.1 


1.82 


1.85 


1.56 


152 


913 


+ 1.7 


3.16 


0.85 


0.28 


107 


513 


+3.4 


5.50 


0.28 


0.10 


44 


174 


+3.2 


9.55 


0.07 


0.04 


17 


51 


+2.8 




— All Data (278.3 


h) — 






1.05 


11.98 


10.80 


283 


1563 


+3.6 


1.82 


1.89 


0.97 


259 


1444 


+3.3 


3.16 


0.64 


0.20 


164 


836 


+3.6 


5.50 


0.29 


0.07 


77 


291 


+4.5 


9.55 


0.08 


0.03 


30 


85 


+3.9 



a: Flux and Errors in units xlO ph cm ^s ^ TeV 
b: Sienificancc from Li&Ma Eq. 17 



using s, b and a normalisation a— 0.143 



(e) Fitted Spectrum: Pure Power-Law 



dN/dE 


= N (E/lTeV)-^ phcm-^s-- 




1999 to 2001 Dataset §1(120.5 h) 


N 


= 4.1 (±2.1stat ± 1.3,y,) X 10" 


7 


= 1.7 (±0.3stat ±0.3syB) 




— 2002 Dataset §2(157.8 h) — 


N 


= 9.3 (±2.9stat ± 1.4sys) X 10" 


7 


= 2.1 (±0.2stat ±0.3,y,) 




— All Data (278.3 h) — 


N 


= 6.2 (±1.5,tat ± 1.3,y,) X 10- 


7 


= 1.9 (±0.1,tat ±0.3sys) 



quasar Cygnus X-3 for which an upper limit {E > 0.7 
TeV) 1.7x10-^"^ ph cm-^s"^ for steady emission is set. 
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Table 3. Summary of point-like flux upper limits for other 
sources in the FoV. Positions for the various X-ray sources 
are taken from Roberts et al. H2U(Jlfl . An energy thresh- 
old of iJth > 0.7 TeV is estimated based on the average 
zenith angle z of observations: i?th = 0.5 cos(z)~^-^. The 
template background model been used [a = 0.168), with 
s. b and S defined as in table □ 



Other Sources: (tight cuts) 9 < 0.12°, w < 1.1, ntei > 3 



Source 


S 


b s 


— ab 


S 


a ^99% 


GeVJ2026+4124 


541 


3429 


-35 


1.3 


0.29 


AXJ2027.6+4116 


771 


4497 


16 


+0.6 


0.42 


GeVJ2035+4214 


784 


4457 


35 


+ 1.2 


0.45 


AXJ2036.0+4218(Srcl) 


620 


4005 


-53 


-1.8 


0.22 


AXJ2035.4-h4222(Src2) 


663 


4083 


-23 


0.7 


0.27 


AXJ2035.9+4229(Src3) 


584 


3626 


-25 


0.9 


0.28 


Cygnus X-3 


857 


5613 


-86 


-2.6 


0.17 


a. <plf° =99% upper limit Eth > 0.7 ToV [xlO" 


Ph 


cm-"s-i] 
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